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Recent intensive studies of light slowdown and storage in optically driven atomic systems have been motivated both 
by their fundamental aspects and potential applications. The minimum model system in which such processes can 
occur includes an atomic medium of three-level atoms irradiated by two laser beams (signal and control one) in the 
A configuration. The underlying physical scenario consists in modifying elcctromagnetically induced transparency 
(EIT) in time and, as a consequence, in a reversible mapping of the laser signal onto the atomic excitation represented 
by a ground state atomic coherence; for recent reviews see, e.g., Refs. [H, A more advanced system is a medium 

of four-level atoms interacting with three laser fields in the tripod configuration. The essential new element is that 
more atomic coherences become relevant and one has more external parameters to control the processes. One of 
the possibilities of exploring such a system is a situation in which two control fields are used during the first stage 
of the storage procedure to map a signal onto a certain superposition of two atomic coherences. The form of this 
superposition is determined by the parameters of the control fields. Due to the possibility of creating the other 
superposition of the coherences, orthogonal to the former one, another pulse can be stored at the following stage, 
without affecting the excitation created at the first stage. Similarly the stored excitation can be released in two 
steps, in general from different orthogonal superpositions of the coherences. This means that our system is capable of 
mixing two signals separated in time and of releasing the signals at two different time instants (signal mixing and/or 
splitting) . One can thus make two time-separated pulses interfere at the storage stage and study statistical properties 
of the retrieved signals. Note that the situation, in which there is only a single stored signal in the classical description, 
in a complete quantum analysis requires treating the other signal as "incoming vacuum" which does influence the 
outgoing signals. The discussed system can thus be considered a kind of a beam splitter or interferometer in the 
time domain, operating on stored light, and the above-mentioned two consecutive stages of storing (releasing) light 
can be considered input (output) ports. Moreover, additional manipulations on atomic coherences are possible at the 
storage stage. An essential difference compared with standard devices is that four light beams (two incoming and two 
outgoing ones) propagate in the same direction but are well separated in time. In contradistinction to usual beam 
splitters here it is possible to smoothly regulate amplitude and phase relations of the outgoing light (as if one could 
change on demand the transmission of the splitter and the phase jump of the reflected light). With such a flexible 
analogue of a beam splitter it is reasonable and interesting to check how statistical properties of the outgoing light 
depend on the parameters of the control fields and/or additional interactions at the storage stage. Note that some 
statistical properties of photon states transformed by a usual beam splitter have been discussed in Ref. [1] 

A medium of atoms in the tripod configuration has already been discussed in the context of light propagation. 
Paspalakis and Knight [5[ studied propagation and slowdown of a weak classical probe pulse and discussed the 
susceptibility of the medium dressed by two control fields. Wang et al. [|| theoretically performed splitting of a single 
photon by storing it in one atomic coherence, by transferring a part of the excitation to the other one by F-STIRAP 
technique and finally by releasing the pulse from the latter. Mazets investigated a simultaneous adiabatic propagation 
of three laser beams of comparable intensities 7] and has analyzed propagation of two quantum signals propagating 
in the tripod medium dressed by a single control field and demonstrated the existence of two modes propagating with 
different velocities Q . Petrosyan and Malakyan Q considered propagation of two orthogonal polarization components 
of a single beam and have studied the effect of nonlinear interaction of the two probes. In our earlier paper j 101 ] we 
have studied storing of a classical light pulse in the tripod medium driven by two control fields. 

An important class of questions are those concerning quantum statistical properties of light pulses, not only because 
of their fundamental physical aspect but also in the context of information storage and retrieving. In Ref. we have 
inve stig ated quantum interference of photons in this context and we have shown how to realize the Hong-Ou-Mandel 
[H, [U interference in such a system. We have also discussed how certain statistical properties of the stored light 
are modified by vacuum fluctuations. The aim of the present paper is to generalize those investigations in the case in 
which the input light states are various non-classical states, the analyzed quantities are not only the photon number 
but also quadrature operators and the homodyne signal. In particular we study the dependence of the statistical light 
properties on the parameters of the control fields. 

The atomic medium considered is composed of four-level atoms with an upper state a and three lowers states b 
(initial), c and d. The quantum signal field 1 couples the states b and a while two classical control fields 2 and 3 
couple the state a respectively with c and d. 

The realization of beam splitting/mixing consists in storing two light pulses in two consecutive time-separated steps 
by switching off at the same rate two sets of the control fields of precisely chosen amplitudes (slowly varying at the 
same rate) and phases (assumed constant); for details see ref. [ll[. Switching the control fields off makes the medium 
opaque and the signals are stored independently in coherent superpositions of the atomic coherences. In the first step 
the Rabi frequencies and phases of the control fields are arbitrary, say, respectively, 0°, ^3, X2 an d X3: with the mixing 
n n° 

angle is defined as qr = arctan ^ = const. In the second step the corresponding control fields are characterized by 

the angles (? — 4>q, \\ + 7r ; xi)- As a consequence, two signals with shapes given by the wave packets fi(z) and fi(z) 
of overlap s = J f*{z)fi{z)dz are independently stored in two orthogonal combinations of the atomic coherences o~b c 
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and <Jbd- Releasing the pulses takes places again in two time-separated steps by switching on at the same rate two 
sets of the control fields characterized, respectively, by arbitrary angles ((f) 1 , x 2 , x\) an d (f — <^ 1 ; X 2 + tTj Xji)- The 
relations between the amplitudes and phases of the control fields at the storing and release stages determine how, 
i.e. with which amplitudes and phases, the stored excitation is divided into two time-separated outgoing signals. An 
additional control procedure is also possible consisting in modifying the atomic coherences at the storage stage by 
some additional interaction. 

We assume that the evolution of our system is adiabatic and can be described perturbatively with respect to the 
signal field, resonance conditions are fulfilled and relaxation processes and noise effects can be neglected. In fact it 
has been proven by Peng et al. [3 that such quantum properties like squeezing and entanglement do not suffer much 
due to dephasing and noise in EIT A medium. Note however that this observation may not be true in general, e.g., for 
the two beams being of comparable intensities or far from the two-photon resonance [la ] . Under such assumptions the 
evolution can be described in terms of joint atom+field excitations called polaritons which are described by photon 
and/or atomic excitation annihilation operators 

X°{z, t) = {e (+) (z, t) cos - Jfc gmO[exp(iy&) cos (f>°a bc (z, t) + exp(<x§) sin *)]}, 



= ^{^{eMiX^)sincf> a bc (z,t)-ex P {ix 3 )cosct> a bd (z,t)}}. (1) 

In the above formula e^ + ^ is the positive frequency part of the signal electric field, d is the dipole moment operator, 
&af3 = \<x >< /9 1 exp(ik a pz) are atomic flip operators, g = \J 2 ^y , ^ 2 = N\d a b\ 2 uj/(2eQh) 1 N is the atom density, 
V is the quantization volume, L is the length of the sample, u> is the field frequency, k a p = (E a — Ep)/(hc), 
tan 2 9 = k 2 /[(SI-]) 2 + (^Is) 2 ]- Changing the control fields corresponding to the first stage of storing (release) to those 
corresponding to the second stage means in fact an interchange of the polaritons X± and X 2 . 

Thus an annihilation operator of a photon stored at the first (second) stage as a packet /i (f 2 ) becomes an atomic 
annihilation operator A°(l) [X 2 (2)), obtained by expanding the operators of Eq. (1) in the two-element basis (in 
our previous paper [ll| we have denoted X\ by ^> and A 2 by Z). In an analogous way the operators of the outgoing 
photons, before they have been released, are atomic excitation operators Xj 2 (1)2). The main idea of the theoretical 
description is based on representing each polariton operator at the release stage by a superposition of operators at 
the storage stage (or vice versa) 

Xj(m) = ]T S jk X° k (m), (2) 

fe=l,2 

where j — 1,2, m — 1,2 and the unitary matrix S 

S u = cose/) 1 cos</>°exp[i(x2 - xl)] +sm</> 1 smcjp exp[i(x£ - x%% 
S 12 = -cos^ 1 sin0°exp[i(x2 - xl)\ +sin<^ cos cfP exp[i(x3 - X°)L 
S 2 x = -sin^ 1 cos0°exp[i(x2 - X°)] +cos</) 1 sin^ exp[ix3 - x°)}, 
S 22 = sin^ 1 sin</>° exp[i(x2 - X 2 )} + costj) 1 cos<f>° exp[ix3 - X°)]- 

(3) 

The commutation relations pertain in general to nonorthogonal modes (due to an overlap of the wave packets fi f2 ) 
and read 

[XUm),Xf(n)]=S jklmn , (4) 



where j, fc, to, n = 1,2 , p = 0,1 and 



7=i in- ( g ) 



This formalism allows one to investigate how quantum statistical properties of the released light are modified due 
to its storage, depending on the characteristics of the control fields and possibly on phase operations on the atomic 
coherences into which the pulse has been mapped. 

Consider first the case of quantum interference of multiple photons [16fl . The incoming photons states are then 
definite Fock states. Note that the cases of interference of three and four photons on standard, possibly non-symmetric 
beam splitters have recently been studied both theoretically and experimentally (l7l . [iH ] . 
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Storing n and m photons at, respectively, first and second storing stage means that the medium is in the quantum 
state 

|4, >= ^=(X 1 t(l))"(X 2 °t(2)) m |0 >, (6) 
V n'.ml 

where the "vacuum" |0 > denotes the state in which all the medium atoms are in the ground state. Expressing the 
operators X® by X]. (Eq. (2)) one gets 

>^ ^ >= E E (fc) [i S h ll S^ k S\^S^ +k [xl\l)] k [xl\l^ > 

Each component corresponds to releasing i photons at the first release stage and n + m — i photons at the second 
one, irrespective from the spatial localization of the excitation. The norm < tpi\ipi >= P(i) gives thus the probability 
of such a photon final distribution, which can be easily calculated when storing at the two stages occurred at the 
same localization and with the same profile fx = f 2 (i.e. s = 1; in such a case we drop the arguments 1 and 2 of X) 

Pii) = ^^E© (£) (.^) (. TO fc/ )^^^l-5™--^ r r5- fc '5«-^™-^'. (8) 

For non-overlapping packets inside the sample, i.e. for s ^ 1, the calculations are tedious and can be performed in a 
recursive way. 

One can also get simple formulae for the mean value and variance of the final photon distribution. The polariton 
number operator to describe photons released at the first release stage in the case of arbitrary packet overlap s is 

N ' = EC-^^OKKfcfri*- (9) 

1 1 Jk 

The latter formula can be obtained as a sum of numbers of excitations from Schmidt-orthogonalized modes and 
reduces to the standard expressions: X\(l)Xi{l) for s = 1 or X\{1)X\(1) + X\(2)Xx(2) for s = 0. An analogous 
formula holds for the photons released at the second stage described by N%- The mean values of and its square 
can be calculated from Eqs. (9), (6) and (2) by a multiple usage of the commutation relations to finally give 

N} = \S 11 \ 2 n+\S 12 \ 2 m, (10) 

and 



[Nff = |5 n | 4 n 2 + \S 12 \ 4 m 2 + \S n \ 2 \S 12 \ 2 [2nm(l + \s\ 2 ) + m + n]. (11) 

The Fano factor being the ratio of the variance W and the corresponding mean value is thus 

W{N\) _ \S 11 \ 2 \S 12 \ 2 (2nm\s\ 2 + m + n) 

|SH| 2 n + |S 12 | 2 m ' 1 j 

In particular for equal photon numbers in the incoming channels (n = m) 

S^^2|5 11 | 2 |5 12 | 2 H S | 2 + 1). (13) 
Nl 

Eq. (12) is a generalization of our earlier results obtained in the case of vacuum in the second input channel. Though 
the variances in both input channels are zero the variance in the output channel 1 (as well as in the output channel 
2) is nonzero. Thus the outgoing states are no longer Fock states and the Fano factor, being here the measure of 
deviation of the ouput state with respect to the input state, is maximum if the overlap s = 1, i.e. when both input 
signals have the same shape and have been stored at the same place inside the sample. One can also see a continuous 
transition to the situation in which the input signals do not overlap at all. 

It can be seen from Eq. (1) that instead of changing the phase X2,3 of the control fields one can realize an equivalent 
modification of the polaritons by introducing additional modifications of the atomic coherences cr a g. For example, as 
we have mentioned in our previous paper on the propagation and storage of classical pulses fTo| , one can make such a 
choice of the atomic states that switching on a magnetic field parallel to the propagation direction induces level shifts 
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of the atomic levels and thus a phase shift 5 of the coherence a\,c while leaving Uhd unchanged. If, e.g., the states 6, 
c, and d are Zeeman sublevels of the hyperfine structure with the quantum numbers (F, M) taken respectively (2,-1), 
(2,1), and (1,1), the coherence <7f, c will acquire a phase shift of absolute value 6 = eA/(2m), where A is the area of the 
magnetic pulse (A — Bt for a rectangular pulse of induction B and duration t). Introducing the phase 8 is equivalent 
to changing the phase x 2 . Thus if the amplitudes of the control fields 2 and 3 both at the storing and release stages 
are equal, which means (ft = (ft 1 = 7r/4 and the phases of the control fields are zero, the polariton transformation 
matrix (Eq. (3)) has the elements Sn = S22 = exp(— i|) cos |, S12 = S21 = i exp(— i|) sin 4. An experiment of the 
Hong-Ou-Mandel type would then consist in trapping single photons (m = n = 1) of the same shape at the same 
place of the sample. The probability of releasing a single photon at the first release stage (i = 1) given by Eq. (8) is 

P(l) = cos 2 <5, (14) 

which means that if no phase modification has been made, the two channels do not mix at all and each photon is 
released as it has been trapped. The effect can be modulated by changing the phase 6 due to a magnetic pulse. Photon 
coalescence occurs for 6 — \. Another typical realization of the standard HOM interference occurs for \\ = 0; in such 
a case P(l) = cos 2 (0 1 — (ft ). 

As a more general example we illustrate in Fig. [1] the situation in which n = m = 6 photons have been trapped 
in each input channel, in the case of the amplitudes of the control fields at the storing stage such that (ft = it/ 8 and 
the corresponding zero phases (x 2 = Xa)- We show the probability of obtaining exactly 6 photons in the first output 
channel for the phase x\ — as a function of the ratio of the amplitudes of the control fields at the release stage 
(represented by (ft 1 ) and the phase of the first control field x\ ■ Apart from 3 peaks reaching unity, which correspond 
to the situations when the two channels do not mix at all or are interchanged, one can see such combinations of the 
parameters (ft 1 and x\ that the probability drops to zero. Such a situation can be called generalized Hong-Ou-Mandel 
interference and consists in zero probability of the photon configuration in the output channels identical with that in 
the input channels. It is however known that a complete photon coalescence in any of the two output channels can 
occur only in the case of storing single photons in both input channels [l3l ] . Increasing the number n of photons in 
the two input channels leads to an increase of the number of local minima (zeros) of P(n). 

Another interesting situation occurs when one stores a quantum squeezed state in each input channel, assuming 
overlapping wave packets (s = 1), which means that after storage has been performed the medium is in the state 
\a\r\ > \0L2r2 > (we consider real squeezing parameters r\^), where 

a n 1 

\a 3 r 3 >=^^ T A]"|0>exp(--|« J | 2 ), (15) 



where Aj = cosh r 3 X° + smhrjXj , The quadratures corresponding to the outgoing light in the first output channel 
are given by 

p=^={Xl-X?). (16) 

Expressing the operators X\ in terms X° 2 and taking advantage of the fact that aj is the eigenvalue of Aj we obtain 
the mean values 

1 / 

q = — p(ui«i + u 2 a2) + c.c. 
v2 

¥ = ^[|wi| 2 (l + 2|ai| 2 ) + |u 2 | 2 (l + 2|a 2 | 2 )] + ^[ujaj + u\a\ + 2uiu 2 aia2 + 2uiu* 2 atxal + c.c], (17) 

where u\ — S'ncoshri — S^sinhfi, and u 2 = Svi cosh r 2 — S 1 ^ sinh r 2 . Analogous expressions hold for p and p 2 
except that S11 and S12 must be replaced by —iSn and — 1S12, respectively. The results are shown in Figs and [3] 
If n = 7'2 > one can easily show that squeezing in q for the outgoing pulses cannot be smaller than exp(— 2r\)/2, 
which is equal to squeezing in q of the input states, and cannot be larger than exp(2ri)/2, which is squeezing in p of 
the input states. One can observe that for real ri j2 both incoming pulses minimize the uncertainty relation which is 
no longer true for the released pulses. The product of variances of the quadratures corresponding to the photons in 
the first output channel is shown in Fig. |4]as a function of the angles (ft 1 and \\. Though the product of variances in 
each input channel is exactly |, the output values essentially depend on the mixing properties of the system. 
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The above-described procedure of storing of two light pulses allows for a realization of a homodyne-like analysis of 
the retrieved pulses due to phase-dependent operations on atomic coherences in which those pulses have been stored. 
Suppose that one stores photons in a squeezed state \a\r\ > (see Eq. (15)) at the first storage stage and a Poissonian 
light in the state |a 2 > at the second stage. This is a case which for a symmetric beam splitter constitutes a basis of 
the balanced homodyne method of detection. The interesting characteristics of the released light is now the difference 
between the number of photons outgoing in the first and second output channels given by the operator 

K = X?X\ - X?Xl = (|5 u | a - ISail 2 )*!**! + (l^i 2 | 2 - |S 22 | 2 )X 2 0t X 2 ° + 

r q* q q* q \ yOf vO i ( q q* q q* \ y^t vO /i q\ 

Note that the variance of the two last components in the r.h.s of Eq. (18), for the field 2 treated classically (i.e. 
with the field operators replaced by c- numbers), is proportional to the variance of the quadrature of the signal field 
1 at the angle determined by the phase of the field 2, the phases of the control fields and/or phases due to additional 
interactions. Thus in a symmetric situation (balanced homodyning) in which \Sij\ = \ which occurs, e.g., in the case 
of zero second control field at the storing stage ((ft — 0), equal amplitudes of the control fields at the release stage 
((f) 1 = j) and zero phases — X§ = x\ = we obtain 

K = exp(-ixl)X?X% + h.c. (19) 

We again express the operators X® and XiOf in terms of Ai and A[ and take advantage of the facts that the eigenvalues 
of Ai and X 2 are, respectively a.\ and a 2 . That variance of K is for a classical field 2 of phase 7 

W(K) = |a 2 | 2 [cosh2n - sinh 2n cos 2(x^ + 7)]. (20) 

The detection of the squeezed state in the first input channel is thus sensible equivalently to the phase of the probe 
classical field in the second input channel (as in the usual formulation for a beam splitter) as well as to the phase of 
any of the control fields. As discussed above, changing the phase of a control field may be equivalent to a phase shift 
due to an additional interaction at the storing stage. 

The results are somewhat more complicated in a non-balanced situation. In Fig. [5] we show the variance W(K) as 
a function of the angle (ft 1 and the phase 7 of the Poissonian field in the case of squeezed vacuum in the input channel 
1 . The Poissonian field has been treated quantum-mechanically. Then 

W[K) = cos 2 (2(/. 1 - 2</>°)(isinh 2 2ri + |a 2 | 2 ) + sin 2 ^ 1 - 2</>°)[|a 2 | 2 (cosh2ri - sinh2ri cos2 7 ) + sinh 2 r x \. (21) 

One can see that for our data the deepest oscillations as a function of 7 occur for such amplitudes of the control fields 
at the release stage that (ft 1 = i.e. when (ft 1 — <p° = ^ , so the wave mixing occurs in a symmetric way. For a given 
value of cf) 1 the plot represents two periods of a sine curve with the amplitude of oscillations and their mean level 
depending on — (jr. For (ft 1 — <jr the two input channels do not mix, so there are no oscillations with changing 7; 
then W(K) is the sum of variances of photon numbers in the two input channels. If the Poissonian field is treated 
classically the component sinh 2 r\ in the last term of Eq. (21) is missing. 

We have analyzed some statistical properties of the released light, after two quantum light pulses have been stored in 
a medium of atoms in the tripod configuration. The properties of the outgoing pulses have been discussed, depending 
on the amplitude and phase relations concerning the control fields. For incoming Fock states we have investigated 
the probability distribution of the number of photons released in the two channels. For the corresponding variances 
we have also shown how the photon statistics depends on the overlap of the stored pulses. For stored squeezed states 
we have shown that the variances of the quadratures for the released pulses may drastically change depending on the 
control fields. We have also given the general formulae for the variance of the operator of the difference of the photon 
numbers in the two output channels and discussed the possibility of a homodyne diagnostics. 
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Figures 




FIG. 1: The probability of observing exactly 6 photons in the first output channel in the case of 6 photons entering in both 
input channels as a function of the angles and \2 for (f>° = ty/8 and X2 = X3 = X3 = 0. 



FIG. 4: As in Fig. Hbut for the product W(q)W(p). 
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